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The Raman intensity can be well described by the famous Albrecht equation that consists of A and B terms. 
It is well known that the contribution from Albrecht’s A term can be neglected without loss of accuracy for 
far off-resonant Raman scattering processes. However, as demonstrated in this study, we have found that this 
widely accepted long-standing assumption fails drastically for totally symmetric vibration modes of molecules 
in general off-resonant Raman scattering. Perturbed first principles calculations for water molecule show that 
strong constructive interference between the A and B terms occurs for the Raman intensity of the symmetric 
0-H stretching mode, which can account for about 40% of the total intensity. Meanwhile, a minor destructive 
interference is found for the angle bending mode. The state to state mapping between the Albrecht’s theory 
and the perturbation theory allows us to verify the accuracy of the widely employed perturbation method for 
the dynamic/resonant Raman intensities. The model calculations rationalized from water molecule with the 
bending mode show that the perturbation method is a good approximation only when the absolute energy 
difference between the first excited state and the incident light is more than five times of the vibrational 
energy in ground state. 


I. INTRODUCTION 

Raman spectroscopy^ is one of the standard vibra¬ 
tional spectroscopic tools that has been extensively ap¬ 
plied in different fields. The intensity of the Raman 
scattering can be generally treated by Albrecht’s Raman 
theory^, in which the polarizability is expanded into two 
terms, i.e. the famous Albrecht’s A and B terms^’^ 


^pq — I^pq T ^pq-) (f) 

where p and q represent Cartesian coordinates. If both 
initial and final electronic states belong to the ground 
state, they can be written as^’^ 

Aj,, ^ PeSe’^Qere^ {v^\v^) (v^\P) 


r^g v'^=0 






oo k,0 0 

r_n ^ e'^ v'^ :e9 ~ ^ 


r^g v''"=0 


-1 oo 0 fe,0 




{v^\Qk\v''){v''\v'') 


-1 OO 0 fe,0 

^ ^ Pe9erqere9 /../ 


vAvA{v^\Qk\B 


r^g v^=0 


oo fe,0 0 

E E {vf\vn{v^\Qk\r), 


( 2 ) 


Electronic mail: said@tlieocliem.ktli.se 


where h is the reduced Planck’s constant, |e^) and je’^) 
are electronic ground and excited states, |u*) and \v^) are 
the initial and final vibrational states of |e^) associated 
with frequency \v^) is the vibrational state of je’^) 
associated with frequency , uj is the frequency of the 
incident light, 0 Je^v^:e 9 v^ is the frequency difference be¬ 
tween and |e^)|u*), p^^r is the transition dipole 

moment between |e^) and \e'^) at equilibrium geometry 
(Qo), Pe 9 e^ is the derivative of p^^r with respect to spe¬ 
cific normal mode Qk- Here we should emphasize that 
all modes could contribute as intermediate states in the 
Raman processes. In other words, \v^) represents all pos¬ 
sible combination of all modes, i.e., \v^) = |u[u 2 • • 
where N is the number of vibrational modes. In some 
cases, the denominator in Eq. 2 would be written in terms 
of vertical excitation energy AE^g, which obeys the rela¬ 
tionship AErg = huje'^:e9 +A. Here A is the reorganization 
energy.^ 

The Albrecht’s theory is commonly used to study the 
resonant Raman scattering,^ although it is a general the¬ 
ory for all Raman processes^. It is well known that, for 
strongly dip ole-allowed transitions, A term is dominant 
in the resonant Raman spectra. It is noted that A term 
does not contribute to non totally symmetric modes.^ 
On the other hand, for weakly dip ole-allowed transitions, 
both the A and B terms could significantly contribute to 
the Raman scattering of totally symmetric modes. In 
these cases, the relative magnitude of them could only 
be determined by quantum chemical calculations.It 
is also noted that the effects of higher-order terms were 
also discussed in recent studies. 

For non-resonant Raman scattering processes, we fo¬ 
cus on the zero-frequency limit, i.e., the uj is neglected. 
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In the traditional treatment of the Albrecht’s theory for 
non-resonant situations, it has been often assumed that 
00ervr:e9v^ ~ ujev.^g? As a rcsult, the denominator in Eq. 2 
is independent of the vibrational state \v^). Therefore, 
the summation of could be calculated prior. Because 
of the completeness of \v^) and orthonormality between 
\v^) and the prior summation will result in the A 
term only responsible for the Raleigh scattering and van¬ 
ishing for non-resonant Raman scattering.^ Nowadays, 
this argument has been widely accepted in the theory of 
non-resonant Raman scattering.^’^’^ 2,13 However, the as¬ 
sumption is apparently too strong. If we only consider 
the fact that AE^g ^ hwyr—hwyi — A, in the framework of 
harmonic approximation and expanding the denominator 
by Taylor series, the component related to an arbitrary 
excited state \e^) in Eq. 2 can be rewritten as 
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FIG. 1. Schematic drawing of the single displaced harmonic 
model used to study the contribution of the A term to the 
Raman intensity. Qo (cuk) and Qr (cc^) are the equilibrium ge¬ 
ometries (vibrational frequencies) of ground state |e^) and ex¬ 
cited state je’^), respectively. All other symbols can be found 
in Eq. 3. 


can be rewritten as 


Here we should emphasize the multimode nature of 
and Eor instance, we have the expression of hv^uj'^ = 
Ylk=i schematic drawings for a single vibra¬ 

tional mode of all definitions could be found in Eig. 1, 
where ojk and are the vibrational frequencies related 
to |e^) and je’^), respectively. The summation A^^ over all 
\e^) would return to Apg. If we consider the Stokes shift 
for the fundamental frequency, i.e. = 0 and 
the first two terms in A^^ vanish due to the orthogonal¬ 
ity of the vibrational wavefunctions in the ground state. 
However, the last term survives and should contribute to 
the non-resonant Raman intensity. 

The finite difference method^^ and coupled-perturbed 
method^^ have also been used to compute the intensity of 
non-resonant Raman spectra by directly differentiating 
the electronic polarizabilities, 

apq = {v^\ae,pqy), (4) 


where^^ 
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Again, the zero-frequency limit is also applied in Eq. 6 
and the summation would become apg. In both cases, 
the Raman cross section can be calculated from the po¬ 
larizabilities in Eqs. 1 and 4.^^ 

In this work, we comprehensively investigate the rela¬ 
tionship between Eqs. 6 and 3. Comparing Eqs. 6 and 3, 
we can immediately find that the Bp^ term corresponds 
to the last term in ap^. Thus, the focus of current study 
is the relationship for A term. We first address this issue 
in the model system and show that the A term indeed 
contributes to final results for non-resonant conditions. 
Then, the A term in time frame is briefly discussed. Ei- 
nally, we take water monomer as a realistic example to 
show the importance of including the A term for final 
results. 
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II. MODEL SYSTEM 

The first hint for the corresponding relationship for 
A term between Eqs. 6 and 3 can be found in the case 
of dAErgjdQk = 0. If we consider that both poten¬ 
tial energy surfaces (PESs) of |e^) and \e^) are harmonic, 
this case is equivalent to the situation that the displace¬ 
ment between the two PESs is 0 (AQ = 0). We can 
immediately obtain that the second term in is 0. On 
the other hand, because of the no shift in PESs, either 
{l\v^) or will be 0 due to the parity symmetry.^ 

As a result, the last term is also 0 for Raman scat¬ 
tering. Hence, one could notice that the last term in 
Apg may correspond to the second term in It is 

well known that the case can be realized in the non to¬ 
tally symmetric vibrations.^ Eor such kind of vibrations, 
there always exists at least one symmetric operator that 
makes the geometries of Qq ± 6Qk to be identity and 
then dAErgjdQk = 0. Thus, the reason of A term has 
no contribution to the non-resonant Raman intensity of 
non totally symmetric vibrations is the orthogonality of 
vibrational wavefunctions between ground and excited 
states rather than the orthogonality of two vibrational 
wavefunctions in the ground state. The former reason is 
exactly the same as that for the selection rule in reso¬ 
nant Raman processes.It is worth to note that this 
conclusion can also be generalized to resonant Raman in¬ 
tensities as well as general PESs (such as double well^^) 
of excited states. 

To further confirm the corresponding relationship be¬ 
tween the last term in A^g and the second term in 
the condition of dAErgjdQk 7 ^ 0 should be enforced. In 
the framework of linear coupling model (LCM)^^ where 
it assumes uuk = both terms can be calculated ana¬ 
lytically and give the same result, be., 

y'2lThjlp%^rq^r^gAQ 

AE^ ' ^ 

When AQ > 0, this result is equal to 

2p%^rq%QghjjJk\fS/AEjig, where S is the Huang-Rhys 
factor^^. We should emphasis that the current treatment 
is consistent with Ting’s work^^ which could be traced 
back to Shorygin’s treatment^^ in 1947. However, 
Ting’s work was restricted to the LCM scheme and the 
vibronic coupling (the B term) has not been included. 
Moreover, Ting’s algorithm was only widely applied in 
the two-state model^^^^^ and intrinsically pre-resonant 
conditions^^’^^’^^. In off-resonance conditions. Ting’s 
A term was usually considered to be very small and 
negligible. 

When going beyond the LCM, with the help of the gen¬ 
eral sum rules for Eranck-Condon integrals,^^ we could 
obtain that the last term in A^g and the second term in 
ajjg both equal to 


We have further performed numerical calculations for 
more general condition of dAErgjdQk 7 ^ 0. The consid¬ 
ered displaced harmonic model is shown in Eig. 1. Eor ex¬ 
ample, in the case of ujk = 1600 cm“^, = 1800 cm“^, 

and AQk = Qr — Qo = a.u., both terms give the same 

value of 0.0055704918 in the unit of 2p%^rq%QghjAE‘jg. 
The analytical and perfect numerical agreement not only 
confirm the equivalence of the two terms but also show 
that, for modes with dAErgjdQk 7 ^ 0, the A term can 
have non-zero contribution to the non-resonant Raman 
intensity. Note that, here the corresponding S is around 
0.4,^^ which is reasonable for realistic molecules. Because 
the je’^) is arbitrary in above discussions, a state to state 
relationship between Albrecht’s theory and perturbation 
method is thus established for non-resonant Raman spec¬ 
trum. It is worth to mention that, in the original work of 
Albrecht, a similar relationship has also been discussed. 
However, it only focused on the case of near-resonance 
conditions (see Eq. 22 in Ref. 2). 

III. TIME-DEPENDENT FRAME 

Time-dependent frame of Raman theory has the ad¬ 
vantage of avoiding the sum-over-states of jv^) in Al¬ 
brecht’s theory.Instead, the dynamics of the 
wave packet given by the initial vibrational wave func¬ 
tion times the corresponding electronic transition dipole 
on the excited states is involved.Mathematically, the 
equivalence of these two representations could be readily 
proven by a simple half-Eourier transform.In time 
frame, the “dynamic” and “static” terms are equivalent 
to the Albrecht’s A and B terms, respectively. Here, the 
“dynamic” term arises from wave packet propagation on 
the Born-Oppenheimer surfaces, while, the “static” term 
comes from the coordinate dependence of the electronic 
transition dipole.In other words, the first and second 
terms in the Taylor expansion of transition dipole, i.e. 
Pe^e9 = P%e 9 E p^r^Qk^ associatc with the “dynamic” 
and “static” terms, respectively. 

In non-resonant conditions, the short time approxima¬ 
tion holds. Thus, in the zero-frequency limit, the time 
average polarizability associated with je’^) could be cal¬ 
culated as^^’^^ 

E = ^ j ( 9 ) 

where the per-factor 2 arises from the zero-frequency 
limit, T = 2'Klij\E^^ — and 

= J dtexg {(pflexpiinrt/h) \4>i). 

( 10 ) 

Here Hr is the vibrational Hamiltonian for je’^), the 
damping factor is omitted, = {v^p e9e^i \4>i) — 

and Eav is the average energy on je’^). If the 
assumption^^’^^ 


( 8 ) 


{(\)f\H,r — Eavl^i) = 0 


( 11 ) 
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FIG. 2. Three vibrational modes of water molecule. 


is adopted, the “dynamic” term in Eq. 10 would be 

(■^^) ( 12 ) 

where AE = Eav — E^g^i. As a result, the “dynamic” 
term (the A term) of the Stokes peaks at fundamen¬ 
tal frequencies also vanishes in the time frame due to 
the orthogonality of the vibrational wavefunctions in the 
ground state. On the other hand, the “static” terms 
survives and its value in the first order of Qk is exactly 
equal to B term in Eq. 3 if we considered AE = AE^g 
(see Appendix B in Ref. 23 for details). This conclusion 
is consistent with that in Ref. 23, where the emphasis for 
far from resonance is solely on the “static” term. 

The reason of the vanishing of “dynamic” term in non¬ 
resonance in the time frame would be easily located as 
the assumption of Eq. 11, which is equivalent to the as¬ 
sumption of 0 Je^v^:e 9 v^ ~ ^e^-.e9 in the energy frame. Both 
assumptions neglect the vibrational contributions in ex¬ 
cited states. To obtain the correct result, in the time 
frame, Eq. 9 should be calculated with the inclusion of 
vibrational contributions in Hr, so that the contribution 
of the “dynamic” term (Albrecht’s A term) will be in¬ 
cluded in the total intensity. 


IV. FIRST PRINCIPLES EXAMPLE: WATER 
MONOMER 

Although the contribution of A term for non-resonant 
Raman intensities has been briefly touched,its mag¬ 
nitude for realistic molecules was still controversial^^’^^. 
Therefore, a first principles example is highly desirable 
to address this issue. To the best of our knowledge, there 
were no previous examples that have investigated the 
magnitudes of both A and B terms at first principles 
level for non-resonant Raman intensities. In the follow¬ 
ing, we take water molecule as a realistic example to fully 
examine the contribution of the A term. Water molecule 
has three vibrational modes,which are ZH-O-H bend¬ 
ing (02), symmetric 0-H stretching (-oi), and asymmet¬ 
ric 0-H stretching modes (03), respectively. The vibra¬ 
tional modes of water are depicted in Eig. 2. All required 
excited states except core-hole excitations were calcu¬ 
lated at time-dependent density functional theory level. 
The finite difference method^^ was used to calculated the 
derivatives in a'L . All calculated a'L were summed up 


TABLE I. Calculated Raman scattering factors for water (in 
the unit of A^/amu) contain only the A term (A), the B 
term (R), or all terms (Tot.), as well as analytical results 
(Anal.). The experimental measurements (Expt.) extracted 
from Refs. 15, 40, and 41 also included for comparison. 


Mode 

Symmetry 

A 

B 

Tot. 

Anal. 

Expt. 

V2 


0.14 

1.54 

1.10 

1.10 

0.9±0.2 

Vl 


11.66 

48.32 

104.94 

104.54 

108±14 

V3 

B 2 

0.00 

25.68 

25.68 

25.67 

19.2±2.1 


over all singly excited states to obtain apq and finally 
Raman scattering factors were calculated by^^’^^ 


Sk=^5al+7^l (13) 


where 


2 i^yy)k 

“ 2 ! ~ i^yy)k] 

+ [i^yy)k ~ (^zz)k\ 
+ [{^zz)k ~ {^xx)k] 


+ 6 


i^xy)]^ + {(^xz)k i^yz)k 


(14) 


Here the subscript “/c” represents the individual vibra¬ 
tional mode. By using Eq. 6, we have calculated the 
scattering factors from three cases: only the A term (5 'a), 
only the B term (5 'b), and both A and B terms (S'xot)- 
With these scattering factors, the corresponding Raman 
intensities, Le., Ib and /xot, are calculated as^’^^ 



h Sk _1_ 

87r^cz>/c 45 1 — exp {—hci>k/kBT) ’ 


(15) 


where cq is the vacuum permittivity, z>in and i>k are the 
wave number of incident light and vibrational mode, re¬ 
spectively, h is the Planck’s constant, c is the speed of 
light, /cb is the Boltzmann constant, and T is the temper¬ 
ature. By definition, the interference contribution is cal¬ 
culated by Jxot — I A — Ib- For comparison, the analytical 
Raman scattering factors at the same density functional 
theory level were also computed. 

All calculated scattering factors as well as experimen¬ 
tal values^^’^^ are listed in Table I. As expected, the 
calculated 5'xot are identical to the analytical results and 
in good agreement with the experimental values^^’^^’^^. 
Here a slight overestimation for vs should be attributed 
to the approximate exchange-correlation functionals, 
since, at the CASSCE(10,10) level with the same basis 
set used here, the calculated value (20.4 A /amu) is in 
the range of experimental observation.^^ Eor the vs mode 











5 


-Total 



FIG. 3. Calculated Raman intensities for water (black solid 
line). The contributions from the A term (red dotted), the 
B term (blue dotted), and the interference (black dotted) 
were also shown. The incident light and temperature were 
set to be 514.5 nm and 400 K followed by the experimental 
conditions,respectively. All spectra have been broadened 
by Lorentzian function with a full with at the half-maximum 
of 10 cm“^. 


that belongs to B 2 symmetry, the A term does not con¬ 
tribute to the total intensity, consistent with the discus¬ 
sion mentioned above. For the Vi and V 2 modes, on the 
other hand, the contribution of the A term can no longer 
be neglected. For instance, the calculated Sb oivi is less 
than half of the analytical value. It is interesting to note 
that Sa can only account for about 10% of the analytical 
result and the summation of contributions from Sa and 
Sb terms can not provide the correct answer neither. On 
the other hand, a simple estimation shows that the value 
from (\/Sa + \/5b)^ (107 A /amu) can reproduce the an¬ 
alytic result nicely. It indicates that a strong constructive 
interference between the A and B terms does exit for the 
vi mode and covers almost 40% of the total intensity. 
A minor destructive interference is identified for the V 2 
mode. In this case, the simple summation of Sa and Sb 
gives a value larger than the final one. Contributions 
from different terms for the Raman intensities are also 
displayed in Fig. 3, which emphasizes the importance of 
the interference for the relative intensity. Notice that the 
first excitation energy of water is around 7 eV, which ob¬ 
viously obeys the condition of uj^r^r.^gyi ^ ooev.^g. But 
even though, the contribution of the A term cannot be 
neglected. Overall, although the water monomer is a 
simple case, it is fully adequate to address our argument. 

We notice that simple extension of Eq. 6 by replac¬ 
ing and with 2[AE^^ ± h{uj + '^E)]^ and 

2[AE^^ ± h{uj + ^F)], respectively,^^ was employed to 
calculate dynamic^^ and resonant Raman^^ intensities. 
Here the damping factor F represents the life time of 
the excited state. In addition, the same extension of the 
second term (the A term, labeled as A^’^ hereafter) in 
Eq. 6 has been used to explain the chemical enhancement 
of surface-enhanced Raman spectroscopy.However, 
the accuracy of such extension is questionable because. 



FIG. 4. Relative error of simply extended A term in Eq. 6 
respect to Eq. 2 at different energy of incident light for the 
model in Fig. 1. Here all parameters are rationalized from the 
first excited state of water molecule with the bending mode 
in LGM. Specifically, AQ = 1.35 a.u., ujk — ^'k — 1600 cm“^, 
and uJe^-.eg = 35 0 ;^. Two F of 100 and 800 cm“^ were used. 
The insert figure shows the incident light around the resonant 
situation. 


near the resonant condition, the Taylor expansion for de¬ 
nominator in Eq. 2 may fail and thus \v'^){v^\ cannot be 
summed prior. 

To verify the accuracy of this simplified extension of 
Eq. 6, we computed the relative error of respect to 
Eq. 2 at different energy of incident light for the model 
system shown in Eig. 1. The result is depicted in Eig. 4. 
Here all parameters were rationalized from the first ex¬ 
cited state of water molecule with the bending mode in 
LCM (see caption of Eig. 4 for details). Two typical F, 
i.e. small one of 100 cm“^ and large one of 800 
were used. Around the resonant situation, different F in¬ 
troduces different behaviour for the relative error. For 
instance, the maximum relative error for small F is much 
larger than that for large F. For the specific case, the 
maximum relative errors are even larger than 1400% and 
130% for small and large F, respectively. In the re¬ 
gion of away from resonant condition, the identical be¬ 
haviour for relative errors with different F is observed in 
Fig. 4. Due to the uncertainty of experimental measure¬ 
ments as well as the limitation of the double harmonic 
approximation,^^ we used 15% as the threshold for the 
calculated According to the model calculations, 

we have found that A^’^ is a good approximation when 
\^e^:e 9 “ Ci;| is larger than 5 cj/c. We should emphasize 
that this value is mainly determined by AQ between two 
PESs. Eor example, when AQ = 10 a.u., the range of 
accurate A^’^ is \uJe^:e9 — oo\ > 10 uJk- Apparently, when 
AQ = 0, the parity symmetry will lead to zero A term 
both in Eq. 2 and Eq. 6. 
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V. CONCLUSION 

We have clearly shown that the contributions from the 
Albrecht’s A term are too large to be negligible for the 
non-resonant Raman intensities of totally symmetric fun¬ 
damental vibrational modes, in contrast with the com¬ 
mon wisdom of the field. The widely employed perturba¬ 
tion method is also found to be too crude for evaluating 
dynamic or resonant Raman spectra. Our findings are 
conceptually important for correctly understanding and 
modeling of Raman scattering processes under different 
conditions. 
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